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TECHNICAL FIELD: 

This invention relates generally to wireless communications systems and methods and, 
more specifically, relates to systems that use multiple radio frequency (RF) carriers, and 
even more specifically, relates to multicarrier communications systems and methods, 
such as Orthogonal Frequency-Division Multiplexing (OFDM) wireless systems and 
methods. 

BACKGROUND: 

OFDM offers the advantages of improved downlink system capacity, coverage and data 
rates for packet data services with high spectral efficiency. This is due at least in part to a 
substantially rectangular spectrum occupancy and the ability to achieve a low-cost 
implementation using the Fast Fourier Transform (FFT). OFDM has been exploited for 
wideband data communications over mobile radio channels, high bit rate digital 
subscriber lines (HDSLs), asymmetric digital subscriber lines (ADSLs), and digital 
broadcasting. Reference in this regard can be made to the following publications: J. A. C. 
Bingham, "Multicarrier modulation for data transmission: An idea whose time has come," 
IEEE Commun. Mag., vol. 28, pp. 5-14, May 1990; J. S. Chow, J. C. Tu, and J. M. 
Cioffi, "A discrete multitone transceiver system for HDSL applications," IEEE J. Select. 
Areas Commun., vol.9, pp. 895-908, Aug. 1991; W. Y. Chen and D. L. Waring, 
"Applicability of ADSL to support video dial tone in the copper loop," IEEE Commun. 
Mag., vol. 32, pp.102-109, May 1994; and to W. Y. Zou and Y. Wu, "COFDM: An 
overview," IEEE Trans. Broadcast., vol. 41, pp. 1-8, March 1995. 

OFDM partitions the entire bandwidth into parallel independent subcarriers to transmit 
parallel data streams. The relatively longer symbol duration provides increased immunity 
to inter-symbol interference (ISI). OFDM has received considerable attention as an air 
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interface for evolution of the Universal Mobile Telecommunications System (UMTS) 
mobile radio systems in the Third Generation Partnership Project (3 GPP) forum. 
Reference in this regard can be made to 3 GPP TR 25.892 vO.l .1 (2003-02), "Technical 
Specification Group Radio Access Network: Feasibility Study for OFDM for UTRAN 
Enhancement (Release 6)'\ 

Channel estimation plays an important role in OFDM systems that is essential to bit and 
power allocation and signal detection. Without perfect knowledge of channel state 
information (CSI), an OFDM system may be inoperative or may incur significant 
performance loss. Pilot-based channel estimations are widely used to estimate the channel 
and equalize the received signal for OFDM systems, where the pilot signals are uniformly 
distributed in the time and/or frequency domain (FD) based on coherent time and 
frequency. Reference in this regard can be made to J. J. Beek, O. Edfors, M. Sandell, S. 
K. Wilson, and P. O. Borjesson, "On Channel Estimation in OFDM Systems," IEEE 
VTC95-Fall, pp.8 15-8 19, Chicago, USA, July 1995, and to O. Edfors, M. Sandell, J. J. 
Beek, S. K. Wilson, and P. O. Borjesson, "OFDM Channel Estimation by Singular Value 
Decomposition," IEEE VTC96-Spring, pp.923-927, Atlanta, USA, April 1996. 

A lowpass filter in a transform domain was proposed for suppressing inter-carrier 
interference (ICI) and additive white Gaussian noise (AWGN) by O.Y. Zhao, and A. 
Huang, "A Novel Channel Estimation Method for OFDM Mobile Communication 
Systems Based on Pilot Signals and Transform-Domain Processing," IEEE 
VTC97-Spring, pp.2089-2093, Phoenix, USA, May 1 997. Second-order interpolation has 
been shown, see M. Hsieh, and C. Wei, "Channel Estimation Techniques Based on Pilot 
Arrangement in OFDM Systems," IEEE Transactions on Broadcasting, Vol. 48, No. 3, 
Sept. 2002, to outperform the linear interpolation with piecewise constant approximation, 
as in J. Rinne, and M. Renfors, "Pilot Spacing in Orthogonal Frequency Division 
Multiplexing Systems on Practical Channels," IEEE Transactions on Consumer 
Electronics, Vol. 42, No. 4, Nov. 1996. Time-domain interpolation is obtained bypassing 
the FD channel estimates into the time domain (TD) through IFFT, zero padding, and 
converting back to the frequency domain through FFT (see Sinem Coleri, M. Ergen, A. 
Puri, and A. Bahai, "Channel Estimation Techniques based on Pilot Arrangement in 
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OFDM Systems/' IEEE Transactions on Broadcasting, Vol. 48, No. 3, Sept. 2002). The 
performance of the time-domain interpolation has been evaluated by R. Steele, "Mobile 
Radio Communications," Pentech Press Limited. London, England, 1992. 

SUMMARY OF THE PREFERRED EMBODIMENTS 

The foregoing and other problems are overcome, and other advantages are realized, in 
accordance with the presently preferred embodiments of these teachings. 

In accordance with this invention an interpolation scheme is provided for channel 
estimation that is based on minimum mean-squared error (MMSE) estimates over 
comb-type pilot signals. The TD interpolated channel estimation suppresses, by the use of 
an advanced noise suppression scheme, the noise jitter that spreads over all or 
substantially all of the bandwidth of interest. The original channel estimates in the FD are 
transformed into the TD by an IFFT function, having a length of the number of pilots per 
OFDM symbol, where in one embodiment a predefined threshold on actual power or 
accumulative power is used to minimize the noise jitter over the bandwidth. In a further 
embodiment the channel estimates in the FD are transformed into the TD and the noise 
jitter is suppressed by preserving the channel estimates at the actual tap delays and setting 
the others to zeroes. In either embodiment the noise suppressed channel estimates are 
then transformed into the FD for FD equalization. Numerical results are used to show that 
the system throughput is significantly improved by advanced noise suppression in 
accordance with this invention, as compared to conventional FD linear interpolation. 

In the further embodiment the delay estimates for frame synchronization are used 
directly, and a method transforms the frequency-domain channel estimates into the time 
domain, preserves the channel impulse response in the time domain at the actual tap 
delays, nulls the remainder of the channel impulse results by zeroing, and transforms the 
noise-suppressed channel estimates back into the frequency domain. 

In accordance with the first embodiment of the invention the weighting factor for the 
power or accumulative power threshold is predefined. Preferably the adaptive weighting 
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factor is a tradeoff between noise suppression and estimation accuracy. A larger 
weighting factor can effectively reduce the noise jitter, but the channel impulse response 
at the actual channel taps may be ignored. Conversely, a smaller weighting factor can 
preserve the channel impulse response at the actual channel taps, but the noise jitter may 
be included as well. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other aspects of these teachings are made more evident in the 
following Detailed Description of the Preferred Embodiments, when read in conjunction 
with the attached Drawing Figures, wherein: 

Fig. 1 is a block diagram of a OFDM transceiver that is suitable for practicing this 
invention; 

Fig. 2. shows a comb-type pilot structure; 

Figs. 3-1 1 are each a graph, where Fig. 3. shows TD channel estimates in Eq. (9) with 
normalized power in an ITU PB3 channel with G=9 dB; Figs. 4A and 4B, collectively 
referred to as Fig. 4, show InfoBER vs. G in a PB3 channel and QPSK modulation; Figs. 
5 A and 5B, collectively referred to as Fig. 5, show InfoBER vs. G in the PB3 channel and 
16QAM; Figs. 6 A and 6B, collectively referred to as Fig. 6, show InfoBER vs. G in a 
VA120 channel and QPSK modulation; Figs. 7A and 7B, collectively referred to as Fig. 
7, show InfoBER vs. G in the VA1 20 channel and 1 6QAM; Figs. 8 A and 8B, collectively 
referred to as Fig. 8, show throughput vs. G in the PB3 channel and QPSK modulation; 
Figs. 9 A and 9B, collectively referred to as Fig. 9, show throughput vs. G in the PB3 
channel and 16QAM; Figs. 10A and 10B, collectively referred to as Fig. 10, show 
throughput vs. G in the VA120 channel and QPSK modulation; and Figs. 1 1 A and 1 IB, 
collectively referred to as Fig. 11, show throughput vs. G in the VA120 channel and 
16QAM, where in each of Figs. 4-1 1 the A designated Figure plots the time-domain 
interpolation function with a power threshold (TD-PT) and the time-domain interpolation 
function with an accumulative power threshold (TD-APT) in accordance with first 
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embodiments of this invention, while the B designated Figure plots the time-domain 
interpolation function that is based on channel tap delays (TD-CTD) in accordance with a 
further embodiment of this invention; 

Fig. 12 shows a Table (Table 2) that illustrates the improvement in dB by the preferred 
embodiments (TD-PT, TD-APT, TD-CTD) versus the conventional FD linear channel 
estimation with InfoBER vs. G; and 

Fig. 1 3 shows a Table (Table 3) that illustrates the system throughput improvement in 
- percentage by the preferred embodiments (TD-PT, TD-APT, TD-CTD) versus the 
conventional FD linear channel estimation with Throughput vs. G (G=10dB). 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The comb-type pilot arrangement, where the pilot signals are uniformly distributed within 
each OFDM block, has been considered and is preferred for use in this invention. The 
comb-type pilot system provides improved resistance to fast fading channels, where the 
pilot spacing is much smaller than the coherence bandwidth of the channel. This 
invention provides there is a channel estimation interpolation scheme that is based the use 
of MMSE estimates over the comb-type pilot signals. 

Fig. 1 is a block diagram of an OFDM transmitter 10 and a OFDM receiver 30 that are 
suitable for practicing this invention. It should be appreciated that the OFDM transmitter 
10 may be embodied in a base station (BS) of a multicarrier communication system, 
while the OFDM receiver 30 may be embodied in a user equipment (UE) or mobile 
station (MS), such as a cellular telephone or other type of mobile user device capable of 
receiving digital voice and/or data signals, such as packet data signals, from the OFDM 
transmitter 10. 

A data source 1 2 provides information bits that are encoded in a channel encoder 14, then 
rate-matched and modulated based on an adaptive modulation and coding ( AMC) set in 
block 16. In block 18 the comb-type known pilot symbols are uniformly inserted in the 
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frequency domain, as is shown in Fig. 2, and in block 20 the OFDM transmitter 10 takes 
an //-point inverse fast Fourier transform (IFFT) of the pilot-added data symbols, such as 

b{n) = IFFT{f?(*)} = Z*(*)exp(y Inkn/N) n = 0,1,2, - - ■ , N - 1 , ( 1 ) 



where B(k) is the data sequence of length N . The output of the IFFT block 20 is 
converted from parallel-to-serial (P/S) in block 22, and in block 24 is inserted by the 
redundancy in the form of a guard interval (GI) of length greater than the maximum delay 
spread, such as 

\b{N + n\ „ = -G,-G + l, • -1 
X[n) = Un\ * = 0,1,2,.,*-1 ' (2) 

where x(n) is the transmitted signals and G is the GI length. Finally, the Gl-added IFFT 
output x(n) is up-converted at the carrier frequency and becomes a transmitted signal 26 
that is sent over a frequency-selective fading channel with AWGN. 

The received signal 28 at the receiver 30 of the user equipment (UE) is given by 
r(0 = A(0<MO+«(0> (3) 



where <g> denotes the convolution operation, h{t)=^a t (t)S{t-r t ) is the channel impulse 



response in the time domain, L is the number of paths, a\(t) is the complex channel 
coefficient at the Ith path, r, is the path delay, j(/)is the delta function and n(t) is the 
additive white Gaussian noise. 

In block 32 the GI is removed, in block 34 serial-to-parallel (S/P) conversion is 
performed, and in block 36 each truncated block is processed by an Appoint FFT as 
follows 

y{n) = r(n + G\ n= 0,1,2,.",* -1 , (4) 
Y{k)=FFT{y(n)}=± N £y{n)exp(-j2xkn/N) k = 0,1,2,-..,* - 1 . (5) 

N n=Q 
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Assuming that the bandwidth of each subcarrier is much less than the channel coherence 
bandwidth, a frequency flat channel model can be assumed at each subcarrier so that only 
a one-tap equalizer 42 is needed for each subcarrier at the receiver 30. For simplicity, the 
received pilot signal can be expressed in vector form such as 
Y p =X p Hp+I P +N p , (6) 

where x p = diag\x p {o\ - ,X p (L p - 1)} is the diagonal matrix with known pilot signals, L p is 
the number of pilot signals per OFDM symbol, I P is the vector of ICI and N P is the vector 
of Gaussian noise in pilot subcarriers. The channel impulse responses at pilot subcarriers 
can be obtained based on, by example, MMSE criteria such as 

Hp=^X p Xp' + ^j X?Y P , (7) 

where the average signal-to-noise ratio SNR = £"{x p }/ a\ , andcr* is the noise variance. 

In Fig. 1 block 38 is the channel estimator and block 40 is the unit that performs an 
interpolation of the channel estimator output, in accordance with this invention. There are 
a plurality of different embodiments of interpolation that can be employed, where one is a 
time-domain interpolation function with a power threshold (TD-PT), another is a 
time-domain interpolation function with an accumulative power threshold (TD- APT), and 
a still further embodiment is a time-domain interpolation function that is based on 
channel tap delays (TD-CTD), all in accordance with this invention. The interpolation 
block 40 may be implemented in circuitry, or in software executed by a data processor, 
such as by a digital signal processor (DSP), or by a combination of circuitry (hardware) 
and software. The output of the interpolation block 40 is provided to the frequency 
equalizer 42. The resulting equalized signal is applied to a demodulator and rate matcher 
block 44 that corresponds to transmitter block 16. The demodulated and rate matched 
signal is then applied to a channel decoder 46 that provides a stream of received bits 48 
that correspond to the bits output from the transmitter data source 12. An optional path 50 
for sending automatic retransmission requests (HARQ) from the receiver 30 to the 
transmitter 10 may be provided to accommodate the presence of errors in the received 
bits 48. 

In conventional frequency-domain linear interpolation (see J. Rinne, and M. Renfors, 
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"Pilot Spacing in Orthogonal Frequency Division Multiplexing Systems on Practical 
Channels," IEEE Transactions on Consumer Electronics, Vol. 42, No. 4, Nov. 1996), and 
more specifically in conventional comb-type pilot based channel estimation methods, two 
consecutive pilot carriers are used to determine the channel response for data subcarriers 
that are located in between the pilots such as 

where S=N/L P is the carrier spacing between two consecutive pilots and [ J is the 
function of round towards minus infinity which rounds the input to the nearest integers 
that are less than or equal to the input element. 

In conventional time-domain interpolation, as in Sinem Coleri, M. Ergen, A. Puri, and A. 
Bahai, "Channel Estimation Techniques based on Pilot Arrangement in OFDM Systems," 
IEEE Transactions on Broadcasting, Vol. 48, No. 3, Sept. 2002, and R. Steele, "Mobile 
Radio Communications," Pentech Press Limited. London, England, 1992, the 
time-domain interpolation is a high-resolution interpolation based on zero-padding and 
IFFT/FFT functions. After obtaining the MMSE channel estimates H p , one first 
transforms it into the time domain by Z^-point IFFT, such as 

//(*)= IFFT XH p {k)cxp(j27lkn/L p ) n = 0,\,2,- i L p -I . (9) 

ir=0 

The IFFT preferably has a length of the number of pilots per OFDM symbol. Then the 
time-domain channel estimate h(n) is interpolated by transforming the L p points into N 
points with a FFT function such as 

H{k)= FFT{h(n)}= — Zh(n)exp(-j2nkn/N) k = 0,1,2, 1 . (10) 

N „=o 

where h{n) = \ h ^ n = °' 1,2 '" yLp ~ 1 is the zero-padded time-domain channel estimate. 

[0 n=L p ,-~,N-\ 

In accordance with a first embodiment of this invention, there is provided by the blocks 
38 and 40 of Fig. 1 a channel estimation time-domain interpolation with a power 
threshold (TD-PT). The FD perfect and actual MMSE channel estimates in the ITU PB3 
channel are transformed into the TD as in Equation (9) and power-normalized, as is 
illustrated in Fig. 3. It is noted that G=Ior/(Ioc+N0), and is the geometry value that 
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applies to the distance between the base station and the user equipment, where lor is the 
total base station transmission power, Ioc is the power of the intercell interference and NO 
is the noise power. 

It can be noticed that the noise is spread over the entire time-domain. In this embodiment 
of the invention a first novel interpolation scheme with a power threshold is provided, 
which suppresses the spread noise jitter in the time domain. The time-domain channel 
estimates with power larger than the predefined power threshold are preserved. 
Otherwise, the channel estimates with power less than the predefined power threshold are 
set to zeros as follows, 

*M=f*M ^N 2>i> <, (ii) 

[0 otherwise 

where | | implies the absolute value, P t =pP a u is the power threshold, p is a weighting 

L '"V i 2 

factor, and P aU = £ is the total power of the time-domain channel estimates. As can 

be seen in Fig. 3, the noise jitter over the entire time-domain can be minimized by setting 
to zero those channel estimates whose power is less than the power threshold. The power 
threshold is defined as p=Q.Q\ 5 2 in Fig. 3. The time-domain channel estimates above the 
power threshold are preserved and transformed into the frequency domain, as in Equation 
(10) for the frequency equalization that occurs in block 42 of Fig. 1. 

If the system is designed to operate at a constant signal-to-noise ratio operation point, a 
constant p can be defined that provides sufficient performance. In a time varying 
environment the power threshold can be adapted most optimally based on the noise and 
interference power. The value of p may be defined to be just above the noise floor level. 
Assuming that the noise and interference is Gaussian distributed with variance a 2 , the 
weighting factor can be defined according to the desired confidence level of selecting 
actual channel taps. For example, if p = 3.84<j 2 or p = 6.66a 2 based on properties of the 
Normal distribution, the confidences would be 95% or 99%, respectively. In practice, 
more convenient values can be used, e.g. factor 4 would correspond to a slightly greater 
confidence than 95%. It can be noticed from Fig. 3 that the adaptive weighting factor is 
the tradeoff between noise suppression and estimate accuracy. A larger weighting factor 
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can effectively reduce the noise jitter, but real channel taps may be ignored, and vice 
versa. 

In an alternative embodiment there is provided a time-domain interpolation function with 
an accumulative power threshold (TD-APT). In this alternative embodiment the power of 
the time-domain channel estimates is ordered as follows 

f.=ORDER||i(»f J /i = 0 > l,2,..,L p -l, (12) 

where ORDER {} is the function to sort elements in ascending order and return an array 
of indices. Then the power of the time-domain channel estimates is accumulated from 
largest to smallest and a threshold on the accumulative power is set for zeroing the 
channel estimates outside of the threshold region, such as 

a(„)=J*(0 V Z\h(i m f <pP all , n=L p -l,.~fl 9 ( 13 ) 
[0 otherwise 

where 90% < p < 95% is the percentage factor which is typically defined to be about one. 
Then the time-domain channel estimates with the accumulative power threshold and zero 
setting can be transformed into the frequency domain as in Equation (10) for frequency 
equalization. 

In a still further embodiment of this invention there is provided a technique for 
time-domain interpolation based on channel tap delays (TD-CTD). The channel tap 
delays are estimated for frame synchronization in OFDM system in order to resist ISI 
with GI removal. The inventors have realized that by using the knowledge of the tap 
delays, the channel estimation can be improved significantly. From Fig. 3 it can be 
noticed that most of the power of the channel impulse responses is allocated at the 
corresponding tap delays in the time domain, while the remainder is basically the noise 
jitter. In this embodiment of the invention the block 40 of Fig. 1 uses time-domain 
interpolation with the knowledge of the CTD to suppress the noise jitter. The channel 
estimates in the FD are first transformed into the TD as in Equation (9). Then, only the 
TD channel estimates at the actual channel tap delays are preserved while the remainder 
are set to zeros as follows: 
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h(n) = 



h(n) if we T\ n = 0,-- 9 L p -1 
0 otherwise 



(14) 



where T = {r l } l = \,--,L is the vector of the tap delays. Time-domain channel estimates 
with CTD-based noise suppression can then be transformed into the frequency domain as 
in Equation (10) for use by the frequency equalization block 42. 

The TD interpolated channel estimation with alternative power thresholds (PT, APT), as 
well as channel tap delays (CTD), disclosed above were evaluated in an OFDM downlink 
system, where the simulation specifications are shown in Table 1 . 



Table 1 . Simulation Specifications 
Systems OFDM Downlink 

Carrier Frequency 2 GHz 

Sampling Frequency 3 .84 MHz 

FFT 512 
Guard Interval 36/37 
Subcarrier spacing 7.5 kHz 

Pilot Signals Comb-type with 10% of total basestation power 

Pilot spacing 16 subcarriers / 120 kHz 

Channel Coding 3GPP 1/3 Turbo 

Modulation QPSK & 16 QAM (Pilot only in QPSK) 

HARQ Maximum 3 retransmissions with 7 TTI delay 

HARQ Feedback Signaling Perfect 

Channel Estimation MMSE with alternative interpolations 

Channel ITU PB3 & VA 1 20 

Figs. 4 through 7 illustrate the alternative channel estimations in ITU PB3 and VA120 
channels with QPSK and 1 6QAM modulation by InfoBER (coded Bit Error Rate) vs. G. 
The advanced channel estimations made possible by the use of this invention significantly 
outperform the conventional schemes, and can approach the performance of the perfect 
scheme with a high G value. With a target InfoBER of 10" 2 , and the conventional FD 
linear scheme as a benchmark, the improvements made possible by the embodiments of 
the channel estimations are summarized in Table 2 of Fig. 12, which shows the 
improvement in dB by the preferred embodiments (TD-PT, TD- APT, TD-CTD) over the 
conventional FD linear channel estimation with InfoBER vs. G. 
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Figs. 8 through 1 1 illustrate the alternative channel estimations in ITU PB3 and VA120 
channel with QPSK and 1 6QAM modulation by system throughput vs. G. With G=l OdB, 
and the conventional FD linear scheme as the benchmark, the system throughput 
improvements made possible by the embodiments of the channel estimations are 
summarized in Table 3 of Fig. 13, which shows the system throughput improvement in 
percentage by the preferred embodiments (TD-PT, TD-APT, TD-CTD) versus the 
conventional FD linear channel estimation with Throughput vs. G (G=10dB). 

By way of conclusion, it is known that the accuracy of the channel estimation has a large 
impact on the system throughput of multicarrier transmissions, and the conventional FD 
linear interpolation does not work well in OFDM DL communication especially for high 
velocity due to the fast Doppler and limited pilot overhead. In this invention, novel TD 
interpolated channel estimations with advanced noise suppression are described for 
multicarrier transmissions. In first embodiments the channel estimates in the FD are 
transformed into the TD and the noise jitter is suppressed by the use of the power 
threshold or the accumulative power threshold techniques. Then the noise suppressed 
channel estimates are transformed into the FD for FD equalization. In a further 
embodiment the channel estimates in the FD are transformed into the TD and the noise 
jitter is suppressed by preserving the channel estimates at the actual tap delays and setting 
the others to zeroes. The noise suppressed channel estimates are then transformed into the 
FD for FD equalization. 

The foregoing description has provided by way of exemplary and non-limiting examples 
a full and informative description of the best method and apparatus presently 
contemplated by the inventors for carrying out the invention. However, various 
modifications and adaptations may become apparent to those skilled in the relevant arts in 
view of the foregoing description, when read in conjunction with the accompanying 
drawings and the appended claims. As but some examples, the use of other similar or 
equivalent circuit and mathematical approaches may be attempted by those skilled in the 
art. Further, while described above in the context of the channel estimation being based 
on minimum mean-squared error (MMSE) estimation over the, preferably, comb-type 
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pilot signals, other estimation schemes, such as least squares estimation, can be also 
employed. However, all such and similar modifications of the teachings of this invention 
will still fall within the scope of this invention. 

Furthermore, some of the features of the present invention could be used to advantage 
without the corresponding use of other features. As such, the foregoing description 
should be considered as merely illustrative of the principles of the present invention, and 
not in limitation thereof. 
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